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摘要 : 超 导 托 克 马 克 聚 变 实验 装置 (EAST) 的 极 向 场 超 导线 圈 发 生 失 超时 ， 超 导 
线圈 电阻 变化 引起 的 微弱 电压 信号 变化 被 强 噪声 淹没 。 针 对 该 问题 ， 运 用 快速 信里 叶 
变化 分 析 线 圈 电 压 信 号 的 时 频 特性 。 根 据 分 析 结 果 ， 提 出 运用 经 验 模 态 分 解 (EMD) 
方法 对 超 导 线 圈 两 端 电 压 信 号 进行 分 析 ， 得 到 若干 固有 模 态 分 量 和 余 项 ， 并 得 知 微弱 
信号 的 能 量 主 要 分 布 在 余 项 中 。 该 方法 能 够 消除 环境 有 影响， 检测 出 电压 信号 的 微弱 变 
化 。 
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体 馈 电 与 优化 Application of EMD Algorithm to Analyze Voltage Signals of 
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Abstract: In experimental advanced superconducting Tokamak (EAST) device, 
during the quench process of a superconducting coil, the voltage variation caused by 
change in coil resistance is so small that is submersed in the background noise. To solve 


the problem, the time-frequency characteristic of the coil voltage signals is analyzed by 


鲍 晓 华 ” 男 1972 年 : ce l ; ; a 
受 ， 博 士 生 导师 ， Fourier decomposition, then the coil voltage signals is decomposed into several intrinsic 


开 究 方向 为 电机 噪声 与 振 | mode functions (IMF) and the remainder by empirical mode decomposition (EMD). The 
动 的 理论 与 控制 技术 。 


result shows that the energy of weak signal is distributed mainly in the remainder and the 
weak variation can be extracted from the noise by EMD. 
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1 引言 


超 导 托 克 马 殉 核 聚变 实验 装置 (Experimental 
Advanced Superconducting Tokamak，EAST) 是 我 
国 自 主 设计 的 全 超 导 大 型 非 圆 截面 托 克 马克 装置 ， 
其 极 向 场 磁体 系统 由 14 个 超 导 线 圈 产 生 ， 为 了 控制 
等 离子 体 击 穿 和 位 形 ， 线 圈 需 要 提供 足够 大 的 交 变 
电流 及 极 快 速 的 磁 通 变化 "“ 。 为 了 保证 装置 的 安全 
运行 ， 需 要 在 超 导 线 圈 发 生 失 超 初 期 及 时 检测 到 由 
于 电阻 变化 引起 的 微弱 电压 信号 变化 。 

失 超 是 超导体 运行 的 重要 现象 ， 极 向 场 线圈 采 
用 CICC 型 导体 制 成 ， 失 超 发 生 时 超 导 线 圈 会 突然 
出 现 电阻 ， 并 逐渐 增 大 。 因 为 电阻 的 变化 会 直接 引 
起 电压 的 变化 ， 所 以 对 电压 的 检测 在 众多 失 超 检测 
方法 中 更 直接 快速 。 极 向 场 超 导 线圈 在 工作 时 处 于 
极为 复杂 多 变 的 电磁 环境 中 ， 线 圈 两 端的 耦合 电压 
非常 强 ， 电 阻 变化 引起 的 电压 变化 会 被 淹没 在 看 合 
噪声 电压 中 ， 因 此 ， 消 除 电压 信号 噪声 ， 准 确 地 检 
测 电压 的 变化 趋势 ， 对 于 失 超 信和 号 的 检测 十 分 重要 。 

近 几 十 年 来 ， 国 内 外 专家 学 者 对 微弱 信号 的 检 
测 问题 进行 了 广泛 的 研究 。 从 传统 的 曲线 拟 合 与 平 
请、 锁定 放 大 、 离 散 量 的 统计 平均 、 相 关 检 测 和 自 
适应 哄 声 抵消 等 技术 ， 到 新 兴 的 随机 共振 、 混 沌 阵 
子 、 经 验 模 态 分 解 和 盲 源 分 离 等 方法 ， 计 算 机 技术 
和 现代 电力 电子 技术 的 发 展 为 实现 这 些 方 法 创造 了 
物质 条 件 ""。 文 献 [7] 介绍 了 曲线 拟 合 原理 ， 文 献 
[8] 比较 了 几 种 不 同 的 平 请 降 噪 方法 ;文献 [9] 构建 
了 锁定 放大 器 并 将 其 应 用 于 多 频率 时 变 信号 的 处 理 ， 
消除 了 低频 噪声 的 影响 ， 提 高 了 输出 信 噪 比 ; 文献 
[10] 详细 讨论 了 取样 积分 法 的 原理 ， 文献 [11] 为 了 
消除 时 域 平均 方法 周期 截断 误差 的 影响 ， 提 出 一 种 
采用 时 域 平均 计算 时 的 最 优 分 段 方 法 ， 相 关 函 数 是 
描述 信号 时 域 相似 度 的 度量 ， 文献 [12] 介绍 了 利用 
相关 函数 的 自 相 关 法 抑制 噪声 ， 提 取 微 弱 周 期 信和 号; 
文献 [13-14] 则 分 别 将 互相 关 检 测 法 用 于 激光 远程 
测 距 中 回 波 信 号 和 齿轮 裂纹 故障 信号 的 分 析 中 ， 有 
效 改善 了 微弱 信号 提取 的 成 功率 和 精度 ， 文 献 [15] 
提出 了 可 变 步 长 自 适应 滤波 算法 ;文献 [16] 针对 广 
义 高 斯 噪声 背景 ， 在 能 量 检测 算法 的 基础 上 提出 了 
改进 的 能 量 检测 算法 ， 仿 真 结 果 证 明了 其 在 低 信 噪 
比 条 件 下 错误 概率 明显 减 小 。 

本 文 利 用 经 验 模 态 分 解 (Empirical Mode Decom- 
position，EMD) 算法 对 加 入 失 超 信号 前 后 超 导 线圈 
两 端的 电压 信号 进行 处 理 ， 得 到 不 同 频率 的 本 征 特 
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征 函数 以 及 余 项 ， 并 通过 希 尔 伯 特 一 黄 变 换 计算 得 
到 各 分 量 的 频谱 及 边际 谱 ， 获 得 信号 的 时 频 分 布 。 
通过 分 析 得 知 微弱 的 电压 变化 能 从 EMD 分 解 的 余 
项 中 体现 出 来 。 


2 经验 模 态 分 解 算法 


2.1 经 验 模 态 算法 原理 

经 验 模 态 分 解 算法 是 由 Norden E. Huang 等 提出 
的 在 内 京 模 态 函 数 (Intrinsic Mode Function，IME ) 
概念 基础 上 ， 将 任意 信号 分 解 成 若干 IME 分 量 的 一 
种 算法 。 每 个 IMF 分 量 都 表征 了 原 信 号 的 局 部 特 
征 ， 可 通过 对 其 解析 信号 相位 求 导 龟 取 瞬时 频率 ， 
与 傅 里 叶 变化 相 比 ， 其 自 适 应 时 频 分 析 的 特点 使 分 
量 的 频率 特性 随时 间 变 化 而 变化 ， 因 此 更 适合 分 析 
非 平稳 、 非 线性 信号 。 但 是 EMD 算法 也 有 其 不 足 
之 处 ， 具 体 表 现 为 以 下 几 点 。 

(1) 端点 效应 现象 。 采 用 样 条 函数 对 极 值 点 进 
行 处 理 时 在 端点 处 会 产生 拟 合 误差 ， 误 差 的 累积 会 
“污染 ”整个 序列 。 

(2) 包 络 线 拟 合 问题 。 三 次 样 条 函数 在 拟 合 数 
据 时 存在 过 冲 和 欠 冲 现象 。 

(3) 模 态 混 全 。 由 于 信号 的 异常 变化 导致 分 解 
后 一 个 模 态 分 量 中 存在 多 个 不 同 的 特征 时 间 尺 度 ， 
使 IMF 分 量 失 去 具体 的 物理 意义 。 

2.2 经 验 模 态 分 解 算 法 步骤 

EMD 算法 的 步骤 如 下 : 

(1) 找 出 原 信 号 的 极 大 值 和 极 小 值 ， 用 三 次 样 
条 插值 函数 拟 合 得 到 上 下 包 络 线 。 

(2) 计算 上 下 包 络 线 的 均值 ， 将 原始 信号 减 去 
该 值得 到 新 的 序列 。 

(3) 判断 该 序列 是 否 符 合 本 征 模 态 函数 的 两 个 
准则 ， 吞 满足 则 可 以 作为 一 个 IMF， 否 则 作为 原始 
信号 ， 重复 前 两 步 ， 直 到 得 到 第 一 个 IMF 分 量 。 

(4) 计算 剩余 信号 ， 将 其 作为 新 的 信号 重复 步 
又 (1)、(2)， 直 到 提取 所 有 的 IMF 分 量 。 

原始 信号 可 表示 为 


xD0= 之 G(O+5() (1) 
式 中 ，c(0) 表示 第 ;个 本 征 模 态 国 数 ，”"(O) 表示 余 项 。 
3 ” 极 向 场 超 导 线圈 电压 信和 号 


3.1 超 导 线 圈 电 压 的 产生 
极 向 场 电源 是 托 克 马克 装置 的 子 系统 之 一 ， 在 


201806.00198v1 


chinaXiv 


等 离子 电流 的 产生 和 控制 中 起 着 至 关 重要 的 作用 。 
它 由 12 套子 电源 系统 组 成 ， 每 套 电 源 由 相应 的 控 
制 、 测 量 、 保 护 和 驱动 等 组 成 ， 通 过 控制 电流 的 变 
化 提供 极 快速 的 磁 通 变化 以 维持 等 离子 体 电流 的 大 
小 和 位 形 。 每 套 电源 系统 的 结构 如 图 1 所 示 。 


控制 
子 系统 | 


整流 测 页 
器 子 子 罗 
系统 系 

统 


失 超 保护 | 
子 系统 
图 1 极 向 场 电源 系统 结构 
Fig.l Block diagram of polodial field power System 


在 每 个 工作 周期 中 ,线圈 中 的 电流 都 会 经 历 上 
升 段 、 平 顶 段 及 对 等 离子 电流 的 反馈 阶段 。 线 圈 两 
端的 电压 、 电 流 由 霍 尔 元 件 测 量 并 传送 至 Windows 
平台 。 图 2 为 一 次 正常 放电 实验 期 间 由 霍 尔 电压 传 
感 器 采集 的 电压 信号 及 霍 尔 电流 传 感 怖 采集 的 电流 信 
号 波形 。 采 样 频 率 为 10kHz， 采 样 点 数 为 230 000 点 。 
图 中 线圈 两 端的 耦合 电压 最 高 达到 300V， 而 由 失 超 
发 生 初 期 的 电压 变化 为 数 伏特 ， 检 测 的 趋势 变化 仅 
占 噪声 电压 的 1% 左右 。 


1000 10000 
0 5000 
这 
二 
至 
-1000 0 
-2000 -5000 
-10 -5 0 5 10 15 20 
时 间 /s 


图 2 线圈 电压 及 电流 信号 


Fig.2 The voltage and current of the coil 


对 于 超导体 线圈 ， 由 于 其 零 电 阻 特性 ， 阁 在 无 
干扰 的 情况 下 通 入 电流 ， 其 两 端的 电压 接近 零 ， 此 
时 能 通过 监测 电压 的 幅 值 快速 判断 超导体 是 否 失 超 。 
但 是 在 极 向 场 电 源 系 统 中 ， 存 在 14 饼 相 互 平行 的 超 
导线 圈 ， 各 线圈 之 间 存 在 磁 通 耦合 ， 等 离子 体 电流 
在 产生 或 破裂 阶段 会 产生 干扰 磁场 ， 装 置 帝 体 中 存 
在 的 感应 袍 流产 生 磁 通 耦 合 干扰 。 因 此 ， 在 装置 运 
行 时 ， 超 导线 圈 两 端 存在 大 量 耦 合 电压 ， 这 些 电压 
对 于 失 超 信号 的 检测 是 极 大 的 干扰 ， 失 超 引 起 的 电 
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压 变化 会 被 冯 没 在 干扰 电压 中 。 
3.2 电压 信号 幅 频 特性 分 析 

为 了 了 解 电压 信号 及 失 超 信号 的 幅 频 特性 ， 利 
用 快速 傅 里 叶 变 换 进行 分 析 ， 结 果 如 图 3 所 示 ， 可 
知 电压 信号 中 有 6 个 幅 值 较 大 的 频率 分 量 ， 对 应 整 
流 桥 输出 电流 中 较 大 的 谐 波 分 量 ， 同 时 存在 很 多 杂 
散 分 量 ， 低 频段 的 信号 尤其 丰富 。 


= > a > [= 全 SS SS > > 

2 2 人 SS SS 0 < [= 《人 二 

[= Le De © [= SS Eo Se Ce be 

-i CN cm 寺 wn Ro Ew Do CN = 
频率 /Hz 


图 3 线圈 电压 幅 频 曲线 
Fig.3 The magnitude-frequency characteristic of 


the voltage signal 


4 ”电压 信号 的 EMD 分 析 


由 傅 里 叶 分 析 知 线圈 两 端 电压 信号 频段 宽 且 有 
丰富 的 低频 成 分 ， 是 一 种 非 平稳 信号 。 失 超 发 生 时 ， 
线圈 电阻 呈 线 性 增加 ， 而 且 失 超 事件 可 能 发 生 在 装 
置 运行 的 任何 时 间 点 。 为 了 模拟 失 超 现象 的 发 生 ， 
选取 5 ~ 7s 时 间 段 的 数据 ， 在 原始 电压 信号 上 加 入 
人 为 构造 的 失 超 信号 ， 即 电阻 值 乘 上 实时 电流 值 ， 
取 发 生 失 超时 刻 为 0 时 刻 ， 得 到 的 失 超 信号 如 图 4 
所 示 。 


| 


: 


0.2 0.4 0.6 0.8 10 12 14 1.6 1.8 2.0 
时 间 /As 
图 4 失 超 电压 信号 


Fig.4 The quench voltage Signal 


失 超 信号 是 非 周期 性 信号 ， 近 似 于 一 个 斜坡 信 
号 ， 其 幅 值 由 0V 逐渐 增加 至 5V， 且 线性 释 加 于 
原 信 号 上 ， 按 照 本 文 第 2 节 的 步骤 对 加 入 失 超 信和 号 
前 后 的 电压 信号 进行 EMD 分析， 提取 不 同 频 率 的 
IMF 分 量 和 趋势 项 。 

图 5 为 加 入 失 超 信号 前 后 的 EMD 分 解 结果 ， 
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时 间 /s 
(a) 正常 电压 信号 的 EMD 分 解 


电压 /V 
已 上 
SS 


.2 0.4 0.6 0.8 1.0 
时 间 /s 
(b) 加 入 失 超 电压 后 信号 的 EMD 分 解 
图 5 电压 信号 的 EMD 分 解 


Fig.3 The extraction signals by EMD 
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其 中 图 5a 显示 原始 信号 被 分 解 为 11 个 IMF 分 量 ， 
低频 分 量 中 存在 轻微 频率 混 炙 效应， 这 是 由 电压 信 
号 中 的 突变 点 造成 的 ， 但 不 影响 趋势 项 的 提取 ;图 
Sb 显示 加 入 失 超 信号 后 ， 电 压 被 分 解 为 12 个 IMF 
分 量 ， 说 明 加 入 信号 的 能 量 并 不 完全 集中 在 趋势 项 
中 ， 而 有 部 分 被 提取 到 了 分 量 中 。 将 原始 信号 减 去 
各 IMF 分量 ,得 到 电压 变化 的 趋势 。 图 6 为 EMD 
分 解 结果 频谱 。 从 分 解 的 最 后 一 层 结 果 看 ， 正 常情 
况 下 ， 截 取 时 段 电压 趋势 范围 在 116 ~ 118V， 而 
加 入 失 超 信号 后 ， 电 压 趋势 值 相应 增加 ， 最 大 达到 
121V。 将 趋势 项 的 变化 与 加 入 的 失 超 信号 进行 对 
比 ， 如 图 7 所 示 ， 用 EMD 分 解 出 来 的 趋势 项 的 变 
化 能 很 好 地 跟踪 全 加 到 原始 电压 信号 中 的 失 超 信号 。 
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图 6 EMD 分 解 结果 频谱 
Fig.6 The spectra of extraction signals by EMD 


针对 EAST 极 向 场 超 导 线 圈 耦 合 噪 声 电压 过 
大 、 瘟 没 失 超 信号 的 问题 ， 分 析 了 电压 信号 的 时 频 
特性 ， 提 出 用 EMD 算法 分 解 电压 信号 ， 从 趋势 项 
中 提取 电压 变化 的 方法 。 通 过 分 析 得 到 如 下 结论 : 

(1) 超 导 线 圈 两 端的 耦合 电压 的 主要 成 分 为 
6 脉 波 ， 同 时 还 含有 丰富 的 低频 干扰 电压 ， 通 过 
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电压 /V 


2 
1 
0 02 04 06 08 10 12 
时 间 /s 
(b) 趋势 项 变化 与 失 超 信 号 对 比 
图 7 电压 信号 趋势 变化 


Fig.7 The variation of voltage signal 


快速 传 里 叶 分 析 可 知 ， 主 要 噪声 
1 000Hz。 

(2) 通过 EMD 分 解 能 提取 不 同 频率 特性 的 电 
压 信 号 ， 从 各 分 量 的 频谱 图 看 ， 分 解 结果 与 原 信和 号 
时 频 分 析 大 致 吻合 。 提 取 的 余 项 在 端点 处 有 偏差 ， 
但 总 体能 较 好 地 体现 加 入 失 超 信号 后 电压 趋势 的 


变化 。 


电压 分 布 在 0 ~ 
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